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Ph. The availability of branching fractions for a large majority of D s decays permits 

the prediction of inclusive branching fractions. This is achieved with the help 
Q\ \ of a modest amount of input from an isospin statistical model applied to non- 

resonant multibody D s decays. A systematic uncertainty in these mostly small 
branching ratios is estimated by comparing predictions of this model with those 
of a model involving quark- ant iquark pair production. The calculated inclusive 
branching fractions can be compared with data (for example, from a large sample 
of D+D*~ + D* + D~ obtained by the CLEO Collaboration) and examined for 



specific final states which can shed light on strong and weak decay mechanisms. 
PACS numbers: 13.20.Fc, 13.25.Ft, 14.40.Lb, 12.40.Ee 



in 

> 

00 . I. INTRODUCTION 

(N 

The mechanisms responsible for decays of hadrons containing heavy quarks are of interest 
both as probes of the strong interactions and as sources of information on the underlying weak 
Q\ '. processes. The identification of signatures of new physics in such processes often relies on a 
firm understanding of long-distance (nonperturbative) effects which can masquerade as such 
new physics. 

^ ■ One incompletely understood process in B meson decays is known as "weak annihilation," 

or WA [lj. One wishes to extract the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vub\ from charmless semileptonic B decays. These constitute only 2% of all semileptonic B 
decays, as iKfe/V^I 2 ~ 1% while phase space favors b utv over b — > civ by a factor of 
2. Strategies thus have evolved to extract the small charmless semileptonic fraction. These 
include, for example, the study of leptons with energies Eg beyond the endpoint for b — > civ. 
The WA process can contaminate the endpoint signal: a B + can turn into a soft 1 = hadronic 
system plus a vector bu which can then annihilate freely into iv. (Helicity arguments greatly 
suppress the annihilation of a pseudoscalar bu into iv.) 

The CLEO [2] and BaBar [3] Collaborations have placed upper limits for WA of a few 
percent of charmless semileptonic b decays, while theoretical estimates [I] lie somewhat lower. 
The WA process is supposed to be of order l/wig, where Q is the heavy quark, so it should 
be more visible in charm decays. A suggestion for probing this process in the semileptonic 
decay — > u>i + V£ was made in Ref. [5]. It was pointed out there that by comparing the 
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decays — > u£ + i>g and -D s — > <p£ + V(. with the corresponding hadronic decays — > c<j7r + 
and -D+ — > 0vr + , one could anticipate a WA contribution to the branching fraction Biui^vg) 
of order 10 -3 , nearly an order of magnitude greater than one would expect from the process 
D s — ► <fi£ + vt taking account of oj-cj) mixing. 

The hadronic decays of D s can shed light on the WA process. The contribution of the 
weak process cs — > ud to the decay D s — > co>7r + is forbidden by G-parity. Charmed particle 
decays to VP final states (V = vector, P = pseudoscalar) dominated by the annihilation 
process cs — > ud do not appear to have a consistent description within flavor SU(3) [6]. In 
particular, D s — > ujti + is expected to be suppressed, whereas it is seen with branching fraction 
B(D S -> to7i + ) = (2.5 ± 0.9) x 10~ 3 , while the allowed mode D s -> is "not seen" [7]. 

The presence of certain D s hadronic decay modes containing an u thus could be regarded 
as evidence for different mechanisms for WA [5]. The decay D s — > co>7r + , forbidden by ordinary 
annihilation, may proceed through pre-radiation of the u, whether via violation of the Okubo- 
Zweig-Iizuka (OZI) rule [8] or rescattering. For instance, the D s can dissociate into two-meson 
states such as D^°K^*> + and D^*' + K^° which rescatter strongly to (cs)u while the virtual 
cs state decays weakly to n + . On the other hand, the decay D s — > u;7r + 7r would be a possible 
signature for ordinary WA caused directly by cs — > ud, where the ud current couples to um+n . 
In this case the cs system must exchange at least two gluons with the final state in order to 
overcome helicity suppression. On the other hand, if cs emits a Q = 1 vector weak current, 
which can couple to uit + , a state ss is left over, which can couple to rj. (The phase space for 
r]' is quite limited.) One would then expect to see D s — > um + rj. 

The study of inclusive u> production in D s decays is thus of interest in shedding light on 
mechanisms of weak decay and their interplay with long-distance (nonperturbative) physics. 
Additionally, the study of rj and particuarly rj' inclusive production is relevant because fits 
based on flavor SU(3) have great difficulty in reproducing the large reported branching frac- 
tions B(D+ -> p + r]) = (13.0 ±2.2)% and B(D+ -> p+r/) = (12.2 ±2.0)% 0, preferring values 
a factor of 2 and 4 less, respectively. 

The CLEO Collaboration has studied the reaction e + e~ — ► (D+Dl~ + D^ + D~) at yfs = 
4.17 GeV (a ~ 1 nb _1 ) accumulating a sample of about 600 nb _1 [9]. The dominant decay 
mode of D*^ is to jDf. It then becomes possible to study a cleanly-identified sample of 
tagged D s decays, and thus to obtain inclusive D s branching ratios. 

In the present paper we undertake a theoretical study of inclusive branching fractions in 

decays, augmenting the extensive list of known processes [7] with estimates for unobserved 
modes based on isospin arguments, particularly those employing a statistical model ^0l[TI] [T2]. 
The statistical model is described briefly in Sec. II. Systematic uncertainties in the predictions 
of this model are obtained by comparison with predictions of an alternative model involving 
quark-aniquark pair production [13]. This model is discussed in Sec. III. We then proceed 
from known D s decay modes to others in Sec. IV, showing that we thereby account for the 
majority of D s decays and presenting a provisional table of all D s branching fractions in Sec. 
V. These are then employed to construct a table of inclusive branching fractions in Sec. VI, 
while Sec. VII concludes. 

II. STATISTICAL MODEL 

In the absence of information on branching fractions for certain nonresonant multibody 
modes, we may use a statistical isospin model to relate them to other known modes. Such 
a model may be constructed by coupling the internal subsystems to isospin amplitudes in 
all possible ways and then assuming the reduced amplitudes are equal in magnitude and 
incoherent in phase [TQl, [TT], [T2] . We illustrate this technique with two cases of total I = I3 = 1, 
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KKir and 3n, using particle orders consistent with those quoted for the isospin Clebsch- 
Gordan coefficients in Ref. [7]. 

A. Example of (KKti) i= j 3=1 

1. (KK)ir: We label the reduced amplitudes Af R by l Kk = 0, 1. Then 

A{K+K-.+) = ±=A kR -\a kR , (1) 

A(K°K°«+) = -±=A™- l ~A K \ (2) 



A(i^TTV) = -^4**, (3) 

so assuming incoherent and equal amplitudes A RR , 

\A(K + K~tt + )\ 2 : \A(K°K°ir + )\ 2 : \A(K + k\°)\ 2 = 3:3:2. (4) 

2. (-kK)K: We label the reduced amplitudes A} K by J wK = 3/2, 1/2. Then 

Afr+K+IT) = ^-Al%, (5) 
A(n+K°K°) = -^=Alf 2 + ^A-f 2 , (6) 
A^K+K°) = -^Alf 2 -±=Al%, (7) 

leading again to the ratios (jlj) if A^% 3 , 2 are equal and incoherent. 

3. (7rif)if: We label the reduced amplitudes A] R by I v g = 3/2, 1/2. Then 

A(n+K-K+) = ~^Al% + ^-Al% , (8) 

V3 
2 

A(n°K°K + ) = ~^Al% - -L^ , (10) 



A(k + K°K°) = , (9) 



leading again to (jlj) if Aj£ 3 i 2 are equal and incoherent. 

B. Example of (37r)j = j 3= i 

The only couplings to consider are (7T7r)7r. For example, choosing a particular order, 

A(n+n+n-) = ^A 2 , (11) 

A(7r + 7rV) = -^A 2 + l -A^ (12) 

so if A2 and Ai are equal and incoherent, 

|A(7r + 7T + 7r-)| 2 : |A(7r + 7r°7r°)| 2 = 3:2. (13) 
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Coupling the pions in a different order one can encounter also an amplitude with I ww = 0, but 
the same result is obtained. 

C. Tables of relative branching fractions 

Results for Cabibbo-favored D s decays are quoted from Ref. [TT| for (KKmr) I=l3=1 final 
states in Table fl] and for (37r)j=j 3= i in Table [III Higher-multiplicity D s decays appear to 
account for a very small fraction of the total [7j. Results for singly-Cabibbo suppressed D s 
decays to K + (wr) can be transcribed from Table I of Ref. [12], which applies to a statisical 
average of / = 1/2 and 1 = 3/2 states for D° — > K + (nir) arising from sudu. Here the 
K + {nil) states arise from sdud, which is related to sudu by isospin reflection. The results 
are shown in Table HTT1 

Table I: Statistical model predictions for charge states in (KKmr)i=i 3 =i. 



n(n + ) + n(n ) 





1 




2 


3 


Q(K) 





— or 







— or 


n = 


1 










1 


1/4 


3/8 








2 


1/10 


9/40 


3/20 


3/10 




3 


1/30 


7/60 


2/15 


4/15 


1/6 



Table II: Statistical model predictions for charge states in (nii) I= j 3= i. 



n(n + ) +n(7r~) 1 3 5 

n = 1 1 - " 

2 1 

3 2/5 3/5 

4 1/5 4/5 

5 3/35 22/35 10/35 

6 1/28 12/28 15/28 



Table III: Statistical model predictions for charge states in K + (nir) arising from singly- 
Cabibbo-suppressed D s decays. A statistical average of contributions from I = 1/2 and 
1 = 3/2 final states has been taken as in Table I of Ref. [T2] . 



n(-K + ) 


+ n(ir~) 


1 


2 3 




n = 1 1/2 


1/2 






2 3/20 


8/20 


9/20 




3 3/45 


9/45 


21/45 12/45 
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III. QUARK-ANTIQUARK PAIR PRODUCTION MODEL 



In order to obtain estimates for systematic theoretical uncertainties in predictions of the 
statistical model we present in this section predictions of an alternative model involving quark- 
antiquark (qq) pair production. 

A. Description of the model 

In this model final nonresonant states involving several pseudoscalar mesons are formed 
by the production of uu and dd quark pairs in every possible way, in association with a 
corresponding four-quark process in hadronic D s decays, and in association with the ud weak 
current in semileptonic r decays. One assumes that amplitudes for a given number of quark 
pairs are equal in magnitudes and incoherent in phase. For completeness, we describe below 
quark diagrams contributing to amplitudes for the three distinct classess studied in Section 
II: 

• Cabibbo-favored D s — > nit and t + — > {nit) + u T . Here one fills in qq pairs (q = u,d) 
between the ud pair produced by the weak current. 

• Cabibbo-favored D s — > KKnn. These decay processes obtain contributions from two 
tree amplitudes which are assumed to add incoherently, a "color-favored" amplitude (T) 
and a "color-suppressed" amplitude (C). In the first amplitude (T), involving the weak 
subprocess c — > (ud)s with a spectator s, one fills in with equal probabilities qq pairs 
between the color-singlet ud pair produced by the weak current and between the final 
s and spectator s. In the second amplitude (C), involving the subprocess c — > (sd)u 
with a spectator s, we fill in with equal probabilities qq pairs between the color-singlet 
sd pair in the weak subprocess and between the us pair. 

• Cabibbo-suppressed D s — > Knir. These decays obtain contributions from both color- 
favored and color-suppressed amplitudes. In the first, based on the weak subprocess 
c — > (ud)d with a spectator s, we fill in with equal probabilities qq pairs between the 
color-singlet ud pair in the weak subprocess and between the final d and spectator s. 
In the color-suppressed amplitude, involving the weak subprocess c — > (dd)u with a 
spectator s, one fills in with equal probabilities qq pairs between the color-singlet dd 
pair in the weak subprocess and between the final u and spectator s. 

Color-favored (T) and color-suppressed (C) tree amplitudes, in processes involving pions and 
a KK pair or pions and a kaon, are independent quantities assumed to add incoherently in 
partial decay rates, r = |T| 2 +|C| 2 . Consider two D s decay processes D s — * f\ and D s — > / 2 
within the same class, / = KKrm or f — Knir for a given n. Denoting corresponding pairs 
of amplitudes by (Ti,Ci) and (T 2 ,C 2 ) and assuming for illustration |Ci|/|C 2 | < |Ti]/|T 2 |, one 
obtains the following lower and upper bounds on ratios of corresponding partial decay rates 
or ratios of branching ratios: 

W < ft = rx < IZil!. (14) 
\c 2 \ 2 b 2 r 2 |T 2 | 2 1 ; 

Ratios |Tx| 2 /|T 2 | 2 and |Ci| 2 /|C 2 | 2 are calculated in the next subsection. They provide lower 
and upper bounds on ratios of branching ratios. These bounds, and ratios of branching frac- 
tions calculated in this model for D s — > mi and t + — > {nit) + i> T , will be compared in Sections 
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IV and V with predictions of the statistical isospin model in order to estimate systematic 
theoretical errors in these predictions. 

B. Tables of relative branching fractions 

Results for relative branching fractions obtained in the qq pair production model are shown 
in Tables IIVI fVT. and IVII The upper and lower parts of the last two tables, denoted T and C 
respectively, correspond to color-favored and color-suppressed amplitudes. 

Table IV: Predictions of the qq pair production model for charge states in (tm)j = i z= \. 



n(n + ) + ra(TT-) 1 3 5 

71= 1 1 

2 1 

3 3/7 4/7 

4 1/5 4/5 

5 5/61 40/61 16/61 



Table V: Predictions of the qq pair production model for charge states in (K Kmr) i=i 3= x, for 
color-favored (T) and color-suppressed (C) amplitudes. 



7l(7T + ) + n(7l ) 





1 


2 




3 


Q(K) 





— or 







— or 


T n = 


1 










1 





1/2 








2 





3/10 


1/5 


1/5 




3 





3/22 


2/11 


2/11 


2/11 


C n = 


1 










1 


1/3 


1/3 








2 


1/9 


2/9 


4/27 


8/27 




3 


1/27 


1/9 


4/27 


8/27 


4/27 



IV. FROM KNOWN D s DECAYS TO OTHER MODES 

A. Hadronic modes from — > t + v t 

As the branching ratio B(D S — ► tv t ) = (6.6 ±0.6)% [7] 114] . and r has important decays to 
hadrons, we must take account of such decays. The known r + branching ratios to hadronic 
final states are summarized in Table IVIH with values in brackets inferred by applying the 
statistical isospin model to £>(r + — > (5n) + u T ) = (7.6 ± 0.5) x 10~ 3 [7J. The second error in 
each of the three values of £>/ corresponding to / + = (57r) + is a systematic error obtained 
by comparing predictions of this model with those of the qq pair production model. Also 
shown are inclusive branching ratios of r + decays to 7r + , 7r°, 7r~, K + , and K°. These are 
multiplied by B(D S — ► tv t ) to obtain contributions to inclusive D s decays. One sees that the 
contributions, particularly to secondary tc + and 7r°, are significant. 
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Table VI: Predictions of the qq pair production model for charge states in K + (mr), for 
color-favored (T) and color-suppressed (C) amplitudes. 



V / 


+ n(7r ) 

1 \ / 





1 


2 


3 


T 


n = 1 





1 








2 





3/5 


2/5 






3 





6/22 


4/11 


4/11 


C 


n = 1 


1 











2 


1/4 


1/4 


1/2 






3 


3/37 


6/37 


20/37 


8/37 



B. Secondary particles from 77, rj', 0, and u; 

The decays of the 77 contribute to secondary pions and photons. In Table IVIIII we note 
these contributions for pions. Secondary photons will not be considered here. 

The decays of the 7/ contribute to a number of inclusive final states. These are summarized 
in Table HXl for the main rj decay modes (those with branching fractions greater than 1%). 

The (f) is produced copiously in D s decays. Its decays to K + K~ and K°K° occur with 
different branching fractions because of differences in the limited phase space. Moreover, it has 
non-negligible branching fraction into 7r + 7r~7r° and 777. In Table 1X1 we collect these branching 
fractions and their implications for inclusive signals. 

Although few D s modes involve u, we consider the main modes of u for completeness in 
Table ED 

C. Treatment of neutral kaons 

Often branching fractions are quoted for decay of a zero-strangeness system into XKs, 
where X has a known strangeness. Thus if S{X) = 1, we infer B{XK) = 2B(XK S ), 
B(XK°) = 0, while if S(X) = -1, we infer B(XK°) = 2B(XK S ), B(XK°) = 0. Similarly 
if a state of known strangeness contains a single Kg we double the branching fraction and 
quote it for the appropriate K° or K . Bracketed states in the table of branching fractions 
will indicate modes for which this has been done. 

D. KKnn final states 

We use the fact that D s — ► KKn is dominated by quasi-two-body 4>tt + , K*°K + , and 
K K* + final states to calculate the relative abundances of three charge states, using a statis- 
tical model for the left-over (presumably non-resonant or broad S-wave) component. We start 
with branching ratios quoted in Ref. [7j: B(D S -> 0tt+) = (4.38 ±0.35)%, B(D S -> K*°K+) = 
(3.9 ± 0.6)% (inferred from the quoted value of (2.6 ± 0.4)% with K*° decaying to K + tt~), 
and B(D S — > K K* + ) = (5.3 ± 1.2)%. We summarize in Table [XTTI the contributions of each 
subsystem to the three KKk charge states, as well as to 27r + 7r°7r~ arising from — > 7r + 7r~7r°. 
Here we have used the branching fractions for in Table [Xj, and have assumed K* decays are 
governed by isospin Clebsch-Gordan coefficients. 

The subtotals represent the sums of quasi-two-body contributions to different charge states. 
As B{D S -> K+K-7T+) = (5.50 ± 0.28)% is measured, we add a nonresonant ("NR") or 
broad S-wave background of (0.74 ± 0.51)% to reach this total. The inferred branching ratios 
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Table VII: Branching fractions Bf = B(r + — > f + v T ) and contributions to inclusive pion 
and kaon production. Modes with branching fractions less than 10 -4 are omitted. Brackets 
denote modes with values obtained from £>(r + — > (57r) + z/ T ) = (7.6 ± 0.5) x 10~ 3 |7J and the 
statistical model coefficients of Table [TTl The last line is calculated with B(D S ) = B(D S — > 



t + u t ) = (6.6 ± 0.6)% [7]. Contributions to inclusive K and K production are found to be 
B(t + -> ~K°Xv T ) = (0.52 ± 0.05)%, B(r+ -> K~Xv T ) = (0.140 ± 0.005)%, or 0.03 and 0.01 
when multiplied by B(D S — > tu t ). Columns 3-7 have been rounded off for convenience. 



/ + 


B f (%) 


B(tt + )(%) 


B(tt u )(%) 


B(tt-)(%) 


B(K + )(%) 


B{K")(%) 


7T+ 


10.91±0.07 


10.91±0.07 














TT+TT 


25.52±0.10 


25.52±0.10 


25.52±0.10 











7T+27T 


9.27±0.12 


9.27±0.12 


18.54±0.24 













n no i n nc 


18.Uo±U.12 


u 


n no i n r\c 

9.U3±U.Uo 


n 

u 


(J 




4 48+n oft 




t: . I -±(J) —1—\J . VJKJ 


A 48-|-n Dfi 


n 


n 


TT + 3ir° 


1 04±0 07 


1 04±0 07 


3 12±0 21 


o 


o 


o 


K+ 


0.695±0.023 











0.70±0.02 





K°ir+ 


0.84±0.04 


0.84±0.04 











0.84±0.04 


K+ir° 


0.428±0.015 





0.43±0.02 





0.43±0.02 





K + 2tt° 


0.063±0.023 





0.13±0.05 





0.06±0.02 







0.287±0.016 


0.29±0.02 





0.29±0.02 


0.29±0.02 





K TT+7r° 


0.39±0.04 


0.39±0.04 


0.39±0.04 








0.39±0.04 


K+TT+1T-1T 


0.081± 0.012 


0.08± 0.01 


0.08± 0.01 


0.08± 0.01 


0.08± 0.01 







0.047±0.021 





0.14±0.06 





0.05±0.02 





K°ir + 2ir° 


0.026±0.024 


0.03±0.02 


0.05±0.05 








0.03±0.02 


K+K° 


0.158±0.016 











0.16±0.02 





K + K°7T 


0.158±0.020 





0.16±0.02 





0.16±0.02 





K°K°7T+ 


0.17±0.04 


0.17±0.04 











0.17±0.04 


K+K--K+ 


0.140±0.005 


0.14±0.01 








0.14±0.01 





K°K°TT+TT 


0.031±0.023 


0.03±0.02 


0.03±0.02 








0.03±0.02 


[3vr+27r-] 


0.217±0.014±0.018 


0.65±0.07 





0.43±0.05 








[2vr + 7r-2vr ] 


0.478±0.031±0.020 


0.96±0.07 


0.96±0.07 


0.48±0.04 








[tt+Att ] 


0.065±0.004±0.003 


0.07 


0.26±0.02 











Totals: 
xB(D s ) 




77.40±0.28 
5.11±0.46 


54.28±0.36 
3.58±0.33 


14.79±0.11 
0.98±0.09 


2.06±0.05 
0.14± 0.01 


1.46±0.08 
O.lOiO.Ol 



Table VIII: Branching fractions in rj decays [7] and contributions to inclusive pions. Branching 
fractions less than 1% are omitted. 



r) mode 


B{%) 


B(7T+)(%) 


B(7T°)(%) 


B(7T-)(%) 


77 


39.31±0.20 











3tt° 


32.56±0.20 





97.68±0.60 





7T + 7r~7r° 


22.73±0.28 


22.73±0.28 


22.73±0.28 


22.73±0.28 


7T + 7T~7 


4.60±0.16 


4.60±0.16 





4.60±0.16 


Total 




27.33±0.32 


120.41±0.66 


27.33±0.32 
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Table IX: Branching fractions in rj decays [7] and contributions to inclusive pions, rj, and uj. 
In calculating pion production we use the branching fractions of Table [VTTT1 for 77 and note that 
B(uj -> TT+Tr-TT ) = (89.2±0.7)%, B(uj -> 7r 7 ) = (8.92±0.24)%, /3(cu -> tt+tt-) = (1.53+°;^)% 
(see Table XI). 



rj mode 




B(ir+)(%) 


B(tt°)(%) 


B(7T-)(%) 


Bfo)(%) 


B(w)(%) 


7T + 7r~77 


44.6±1.4 


56.8±1.8 


53.7±1.7 


56.8±1.8 


44.6±1.4 







29.4±0.9 


29.4±0.9 





29.4±0.9 








7r°7r°?7 


20.7±1.2 


5.66±0.33 


66.3±3.8 


5.66±0.33 


20.7±1.2 







3.02±0.31 


2.74±0.28 


2.96±0.30 


2.74±0.28 





3.02±0.31 


Total 




94.6±2.1 


123.0±4.2 


94.6±2.1 


65.3±1.8 


3.02±0.31 



Table X: Branching fractions of </> [7j and their implications for inclusive pion production. In 
addition the branching fractions for kaon and 77 production in decay are B(K + ) = B(K~) = 
(49.2 ± 0.6)%, B{K°) = B(K°) = (34.2 ± 0.5)%, and B{rj) = (1.30 ± 0.03)%. 



cf) mode 


B{%) 


£(tt+)(%) 






K + K~ 


49.2±0.6 











K°K° 


34.0±0.5 













15.25±0.35 


15.25±0.35 


15.25±0.35 


15.25±0.35 


VI 


1.304±0.025 


0.36± 0.01 


1.57±0.03 


0.36± 0.01 


Total 




15.61±0.35 


16.82±0.35 


15.61±0.35 



Table XI: Branching fractions in uj decays [7] and contributions to inclusive pions. 



uj mode 


B(%) 


B(tT+)(%) 


/3(vr u )(%) 


B(n-)(%) 


7T + 7r _ 7T 


89.2±0.7 


89.2±0.7 


89.2±0.7 


89.2±0.7 


7T°7 


8.92±0.24 





8.92±0.24 





7T + 7T~ 


1 crq+0.11 


1 - d °-0.13 





1 - do -0.13 


Total 




90.7±0.7 


98.1±0.7 


90.7±0.7 



Table XII: Contributions of quasi-two-body D s decays to KKir and 2tt + tt tt° charge states, 
in percent branching ratio from D s . 



Mode K + K~ 



71 



K + TTn° 



27T+71-TT 



(j)TT + 

K*°K + 
K°K* + 
Subtotal 
"NR" 
Total 



2.15±0.17 
2.6±0.4 



4.76±0.43 
0.74±0.51 
5.50±0.28 



1.49±0.12 


3.53±0.80 
5.02±0.81 
0.74±0.51 
5.76±0.96 





1.3±0.2 
1.77±0.40 
3.07±0.45 
0.49±0.34i^| 
3.56±0.67 " 



0.67±0.05 



0.67±0.05 


0.67±0.05 
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Table XIII: Contributions of quasi-two-body D s decays to KK2tt charge states, in percent 
branching ratio from D s . Not shown: A contribution of D s — > <pp + to the 2tt + 7t~2it final state 



with branching fraction (4.0 



■ (15.25 ± 0.35)%/(49.2 ± 0.6)% = (1.24 ± 0.36)%. Totals 



with smaller errors are observed branching fractions [7j and are used instead of model esti- 
mates. (Branching fractions to K + K° 7r + 7T~ and K°K~ 2tt + are inferred by doubling reported 
branching fractions involving Kg-) 



Mode K + K~ 



7T 1 7T 







K"K 7T + 7T 1 ' 



K°K-n + n + K+K 



4.0+li 
1.56±0.56 
5.56±1.28 
0.04±1.37 

5.6±0.5 



TV TV 



Subtotal 
"NR" 
Total 



2.76±0.79 
1.56±0.56 
4.32±0.97 
0.04±1.37 
4.36±1.68 







0.05±1.83±0.02 
1.92±0.26 





3.11±1.11 
3.11±1.11 
0.03±0.91 
3.28±0.24 





0.78±0.28 
0.78±0.28 
0.02±0.61 
0.80±0.67 



B(D S -> K°K°7c + ) = (5.76±0.96)% and 8{D S K + K°ti ) = (3.56±0.67)% (where the error 
includes a small contribution from a systematic model uncertainty) are not very different from 
the respective predictions of (5.50 ± 0.28)% and (3.67 ± 0.19)% of a pure statistical model. 

Observed quasi-two-body decays of D s contributing to KK2n final states are to <p(— > 
K + K~)p + , with branching fraction (4. 01^2)%, and K* + K*°, with branching fraction (7.0 ± 
2.5)%. As seen in Table [XTTTl these two contributions almost completely saturate the final 
state K + K~7c + tt , leaving a very small "NR" amount to be treated model-dependently. 

The only subsystem for which one has reliable information in D s — > KK3tc is 027r + 7r~, 
with quoted branching fraction B(D S — > (f>(— > K + K~)2ti + it~) = (0.59 ± 0.11)%. Using the (j) 
branching fractions in Table [X], one infers 

B(D S -f <p2n + ir-) = (1.20 ± 0.22)% , 
B(D S -> (/>(-> K + K-)2tt + tv-) = (0.59 ± 0.11)% , 



B{D S -> 0(-> K°TT)2n + n-) = (0.41 ± 0.08) 



B(D s ^(f)(^7i + 7i-7i )27i + 7i-) = (0.18 ± 0.03)% . (15) 

The 7r + 7r°7T system accompanies the 2/3 and 3/4 as frequently as 2tt + tt~ in the statistical 
model and the qq production model, respectively (see Tables HT1 and HVT) . Thus one gets the 
corresponding predictions 

B(D S -> 07T+27T ) = (0.80 ± 0.15 ± 0.10)% , 
B(D S -f K + K-)n + 27r°) = (0.39 ± 0.07 ± 0.05)% , 
B(D S ^ K K°)n + 2n ) = (0.27 ± 0.05 ± 0.04)% , 
B(D S -> 0(-> 7T + 7r"7r )7r + 27r°) = (0.12 ± 0.02 ± 0.02)% , (16) 

where the second error represents the difference between predictions of the two models. The 
only D s — > KK37T modes for which branching fractions have been measured are — > 
K+K-2TT+H-, with B = (0.88 ±0.16)%, and D s -> K 8 K 3 2n + n-, with i3 = (8.4 ±3.5) x 10" 4 . 
It may be risky to use the latter to extrapolate to the K°K°2tt + tt~ mode if the kaon pair 
favors one CP eigenstate over another. Consequently, we shall use the difference between 
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Table XIV: Contributions in percent of 4>(3tt) and model contributions to charge states in 
D s — > KK3n . The latter contributions, based on the statistical model, contain a statistical 
error and systematic error obtained by comparison with the qq production model. 



State 



0(3 



7T 



Model 



Total 



K + K-2tt + tt- 
K°K°2tt + it- 

K + K°TT + n-TT° 
K°K-2n+n° 
K + K-7i + 2n° 
K°K°n + 2n° 
K + K°3n° 



0.59±0.11 
0.41±0.08 




0.39±0.09 
0.27±0.06 




0.29±0.19 
0.29±0.19 

0.46±0.30±0.17 
0.23±0.15±0.06 
0.20±0.13±0.02 
0.20±0.13±0.02 
0.06±0.04±0.06 



0.88±0.16 
0.70±0.21 
0.46±0.34 
0.23±0.16 
0.59±0.16 
0.47±0.15 
0.06±0.06 



B(D S — > K + K~2tt + 7t~) and its contribution from the to estimate the amount of the KK3tt 
final state to be treated model-dependently. The results are shown in Table IXIVl where 
symmetric systematic errors are quoted using the maximum of positive and negative errors. 

E. Multiple pions from <p(nn) final states 

We have already calculated branching fractions for D s to decay to multipion states via the 
mode — > 7T + 7r~7r°, but summarize the results here. 

(1) The state 27r + 7r~7r° receives a contribution from D s — > 07r + , resulting in B(Df — > 
2tt + tt-tt°) = (4.38 ± 0.35)% • (15.25 ± 0.35)% = (0.67 ± 0.06)%. 

(2) The state 27r + 7r~27r° receives a contribution from D s — > <pp + , resulting in B(Df — > 
2tt + tt-2tt°) = (4.0lJi)% • (15.25 ± 0.35)%/(49.2 ± 0.6)%, resulting in B(D+ -> 2ti + tt-2tt ) = 
(1.24 ±0.36)%. 

(3) The state 37r + 27r _ 7r° receives a contribution from D s — ► 027r + 7r _ , resulting in B(Df — > 
3tt + 2tt-tt°) = (0.18 ± 0.03)%. 

(4) The state 2tt + 7i~3ii receives a contribution from D s — > 07r + 27r°, resulting in B(D+ — > 
2n+n-3n°) = (0.12 ± 0.03)%. 

F. Avoiding double counting from D s — > (rj, r}') + (ir + , p + ) 

The decays D s — > (7]7i + ,r]'7i + ,rip + ,ri'p + ) can populate numerous multi-pion final states. 
Of these, the only inclusive branching fraction explicitly quoted in Ref. [7] is B(D S — > 
37T + 27r~7r°) = (4.9 ± 3.2)%. We estimate using quoted branching fractions that B(D S — > 
ri'n + )B(r]' -> rin + it-)B(ri -> tt+tt-tt ) = (0.39 ± 0.04)%, leaving (4.5 ± 3.2)% to be accounted 
for elsewhere. It is this figure that we include in the table of D s branching fractions. 

V. PROVISIONAL TABLES OF BRANCHING FRACTIONS 

We summarize D s branching fractions to leptonic and semileptonic final states in Table 
IXVI and to hadronic final states in Table IXVI1 Averages in the literature j7] have been 
supplemented where necessary with statistical model estimates, but an attempt has been 
made to use as much input from data as possible. We treat 7]X, rj'X, and uX final states 
explicitly, but quote only final states for <pX. The sum of the entries in Tables IXVI and IXVII 
is (102.7 ± 5.3)%. The largest source of error is a rather early measurement [16], B(D S — > 
37T + 27r~7r°) = (4.9±3.2)%, which includes a small contribution from D s — > r/7r + . Replacement 
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Table XV: D s branching fractions to leptonic and semileptonic modes [TJ [HI [15] . Only modes 
contributing to inclusive hadron production are shown. Modes with B < 10~ 3 are omitted. 
Each D s — > Xl + vn represents the assumed sum of Xe + v and Xp^v modes and is doubled 
from the value quoted in Ref. [7]. 



Mode / 


B{D S -> /) (%) 


Remarks 




6.6±0.6 


See Table VII for r + hadronic decays 




5.8±1.2 


See Table IVIIII for 77 hadronic decavs 


r}'E+v e 


2.04±0.66 


See Table [IX] for 7/ hadronic decays 


<j>Pv t 


4.72±0.52 


See Table [X] for hadronic decays 


Total 


19.16±1.58 





of this number and many estimated numbers with observed ones will be possible given the 
excellent particle identification and electromagnetic calorimetry of the CLEO-c detector. We 
do not extrapolate to other 6n modes from 37r + 27r~7r°, preferring to wait until it is better- 
measured. 

The last 11 entries in Table IXVT1 sum to a branching fraction of (4.41 ±0.78)% for hadronic 
Cabibbo-suppressed decays of D s , consistent with a fraction \V c d/V cs \ 2 of the Cabibbo-favored 
hadronic decays which amounts to (4.24 ± 0.27)%. 

VI. PREDICTED INCLUSIVE BRANCHING FRACTIONS 

With the branching fractions in Tables IXVI and IXVI| it now becomes possible to calculate 
inclusive particle yields. These are summarized for pions in Table IXVII| for kaons in Table 
IXVIII| and for 77, 77', 0, and u in Table IXIX1 Errors in tt + and tt~ inclusive branching 
ratios can be reduced considerably by improving the branching fraction measurement for 
D s -> 3n + 2TT-Ti . 

VII. CONCLUSIONS 

We have calculated the inclusive branching fractions of D s mesons to several species, 
using the fact that the observed branching fractions, together with modest assumptions about 
unseen charge states, account for all the D s decays to an accuracy of about 5%. Calculations 
of branching ratios for unseen modes, based mostly on a statistical isospin model, involve 
small systematic theoretical uncertainties estimated by comparison with a model using quark- 
antiquark pair production. While many aspects of this analysis bear some resemblance to an 
itemized tax return, several notable features have emerged. 

• The greatest errors on extracting inclusive branching fractions from exclusive modes are 
due to a few final states, notably 37r + 27r~7r°, i]p + , and rj'p + . Improvement of information 
on these modes would be very helpful. 

• The large predicted values for the inclusive branching ratios B(r]) and B(r]') may be 
helpful in determining whether, as suspected in SU(3) fits (see, e.g., Ref. [6j) and in 
factorization calculations or modifications (T7J [TS] based on the observed semileptonic 
decays D s — ► rftiv, the branching fractions for D s — ► r/p + and particularly r]'p + are too 
high. It could be possible that part of the large "signals" for D s — > 7]^p + come from 
misidentified kinematic reflections from other final states. 
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Table XVI: D s branching fractions to hadronic modes [7J. Values for bracketed modes are 
inferred from arguments in Sections II and III. Second errors are systematic uncertainties 
obtained by differences between predictions of the statistical isospin model and the model 
based on qq production. Modes with B < 10~ 3 are omitted. 



Mode / 


B(D S -> /) (%) 


Remarks 


, i), 

[K+Tt] 


2.98±0.18 


Doubled from quoted K + K s value 


K+K-71+ 


5.50±0.28 




_ 

\k°k\+} 


5.76±0.96 


Quasi-2-body + model for NR 


[K + K\°] 


3.56±0.67 


n n 


K + K-tt + it 


5.6±0.5 




[K°K°n + n°] 


4.36±1.68 


Quasi-2-body + model for NR 


[K + K°it + it-] 


1.92±0.26 


Doubled from quoted value with K s 


[K°K-7T + -K + ] 


3.28±0.24 




[K + K\°7l°] 


0.80±0.67 


Quasi-2-body + model for NR 


K + K~7T + 7T + 7T~ 


0.88±0.16 




[K K°n+n+n-] 


0.70±0.21 


0(37r) + model 


[K + K°n + n-n°] 


0.46±0.34 




[K°K~ii + tt + it ] 


0.23±0.16 




\K + K-n + TT°ir°] 

L J 


0.59±0.16 




\K K°7r + ir iT } 


47±0 15 


» » 


\ T<r+W°TT° <7r°<7r°l 
\J\ J\ 7T 7T 7T I 


u.uozcu.uo 


n n 




i 1 1 _i_n ns 

1 . 1 IlCU.Uo 




r„+„o_oi 

I7T 7T 7T I 


u. 1 4tzcu.uozcu.uy 


iviouei 11 om 71 vr 71 


I7T 7T 7T 7T I 


U.U ( ICU.UU 


nuiii <^7T 


r 7 7r+ 


I.OoICU.ZI 






n oK-i-ct no 
u.zozcu.uy 




OVT Z7T 


n 8n-i- n no. 
u.oUzc u.uy 




r97r+7r _ 97r°l 

/i n . / 1 




T^ty^t ) 1 /A n - ^ - q n r\ rnnrl £>1 f tym n ^ •7r~^ _i ?'7r 

J. I vjlll U'A-' CXiLiXX HlvJLlCl 11 vjlll O /I ^ /I 


r +4 Oi 

I7T ^t/T I 


U.Z^tZCU.UOICU.UI 


iviouei 11 om ovr zvr 


?7P + 


1 fi-l-9 9 




[07T Z7T 7T J 


4.0ZC0.Z 


ouDiracimg u.oy/o ior tj it 


Z7T 7T o7T 


U.IZzCU.Uo 


From D s — > (^ti^ti^ti® 


'/ 


Q C-J-f) 4 




v'p + 


1 9 9J-9 n 




[X\ 7T J 


U.ZOlCU.UO 


uouoieQ 11 om jv^vt 


j\ 71 


n ns-i-n no 

U.U01CU.UZ 


(Statistical: 0.25i0.03, qq\ unrestricted) 


K + n 


f( 141+0 D31 

W • .It: J l_ W . WO X 






0.16±0.05 




K + 1T + 1T~ 


0.69±0.05 






0.61±0.04ig;i 


Model from K+n+n^ 


[K+n°n°] 


0.23±0.02lg;i 


n n 


[K°7V + 7l + 7l-] 


0.60±0.22 


Doubled from i^s7r + 7r + 7r~ 


[^ + 7r + 7r-7r ] 


1.05±0.39±0.45 


Model from K 0/ k + 7t + 7t~ 


[K°7r+7r 7r°] 


0.45±0.17 


71 >> 


[K + 3tt°] 


0.15±0.06^ 




Total 


83.57±5.05 
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Table XVII: Inclusive yields of pions from various final states in D s decays. 



Mode 


B(%) 








T + V T 


6.6 ±0.6 


5.11±0.46 


3.58±0.33 


0.98±0.09 




5.8 ±1.2 


1.59±0.33 


6.98±1.45 


1.59±0.33 




2.04±0.66 


1.93±0.63 


2.51±0.82 


1.93±0.63 




4.72±0.52 


0.74±0.08 


0.79±0.09 


0.74±0.08 


K + K-it + 


5.50±0.28 


5.50±0.28 








K°K%+ 


5.76±0.96 


5.76±0.96 








K+K°n° 


3.56±0.67 





3.56±0.67 





K+K-n+iT 


5.6 ±0.5 


5.6 ±0.5 


5.6 ±0.5 





K°K°n + n° 


4.36±1.68 


4.36±1.68 


4.36±1.68 





K+K°it + it- 


1.92±0.26 


1.92±0.26 





1.92±0.26 


K°K-it + it + 


3.28±0.24 


6.56±0.48 








K + K°ir°n° 


0.80±0.67 





1.60±1.34 





K + K~7T + 7T + 7T~ 


0.88±0.16 


1.76±0.32 





0.88±0.16 


K°K°7T + 1T + 7T- 


0.70±0.21 


1.40±0.42 





0.70±0.21 




0.46±0.34 


0.46±0.34 


0.46±0.34 


0.46±0.34 




0.23±0.16 


0.46±0.32 


0.23±0.16 





K + K-ir + 7r°n° 


0.59±0.16 


0.59±0.16 


1.18±0.32 





K°K°n + n n 


0.47±0.15 


0.47±0.15 


0.94±0.30 





XV TV Ji Ji /( 


u.uuznu.uu 


n 

u 


P| 1 fi-UP| 1 s 

u. lonzu. lo 


Pi 
u 


•77" "I" TT TT 

ii ii ii 


i 1 1 _|_n no 
1.11 nzu .uo 




PI 

u 


1 1 1 _i_pi pis 

1 . 1 1_LU . uo 


7r + 7r o 7r o 

ii ii ii 


n 7A-\-r\ 1 n 

U. ( ^ZLLU. 1U 


n 74+n 1 n 


1 zl8-|-n 9D 


PI 
u 


TT TT TT TT^ 
ii ii ii ii 


U.U 1 xu.uu 


i QA-i-D 1 9 

I.OtiZIZU. 1Z 


pi fi7-|_pi pifi 
u.u i znu.uu 


pi fi7-|_pi pifi 
u.u 1 znu.uu 


'I 71 


i c;s-|-n 91 
i . ooznu .z i 


9 ni +n 97 


i qpi+n 9^ 


u .^loznu . uu 




U.ZOILU.Uc/ 


D 48-1-0 1 7 

U.^OXU. 1 i 


pi 9^4-0 DQ 

U.ZOHZU.Uc/ 


pi 994-pi no 


Oii Zf/I 


u.ouzn u.uy 


9 40-1-0 97 


PI 

u 


1 «pi-l-n 1 R 

l.UUXU. 10 


£ ii ii 


q nn+n 48 


r nn_|_n Qfi 


u.uuznu.yu 


q nn+n 4s 


ii L ±i\ 


pi 9zi-i-n 


U.i^XU.UO 


pi Ofi-upi i 9 
u. yuxu. i z 


P| 

u 


VP + 


io.uhzz.z 


i « c. k_i_o on 

1U.003ZZ.OU 


90 fiK-L/i fit; 
zo.uoin z ±.oo 


Q C.C.-UPI fiPl 
O.OOJZU.UU 


OTT Z7T 7T 


/l c.-1-Q 9 


lo.ozcy .u 


/[ Ki-LQ 9 

^.onzo.z 


q n-i-fi 4 
y.unzu.^i 


Z/l /i O/l 


u. iznzu.uo 


u.z^nzu.uu 


P| Qfi-|-P| piQ 
U.OUILU.Uc/ 


PI 1 o_i_n 
u. iznzu.uo 


' / " 


q Q_i_n a 


7 QO-i-n 7Q 


^.u 1 znu .Hlu 


q Kn+n 00 


V P 


i o 9_i_o n 




97 91 +4 4fi 

Z 1 .Zl IL^t .^lU 


1 1 £.4 4-1 8Q 
1 1 . o^tnz i.oy 




U.ZOILU.UO 


u.zonzu.uo 


PI 

u 


PI 

u 




n fiR+D D9 


n 

u 


pi piR+n ri9 

u.uoznu.uz, 


PI 

u 




n i 41 _i_n n^l 

U . It: iznu .UOl 


u. u^iznu .u i 


pi I 7_i_pi pi4 

U. 1 1 XU.U 1 ! 


pi n4+n ni 

u.u^iznu.ui 


K+n 1 


u. iuilu.uo 


pi i c.-i-n 

U. lOILU.UO 


pi 9pi-i-n nfi 

u.zunzu.uu 


pi 1 C. 1 p| p|C 

u. ioilu.uo 




u.uyznu.uo 


P| fiQ-|-P| Pl^ 

u.uynzu.uo 


PI 
u 


P| fiQ-l-PI pi 1 ^ 
U.Uc/ILU.UO 




61±0 35 


61±0 35 


61±0 35 





i\ TT TT 


n oQ_i_n 1 s 
U.zoztU.lo 


U 


U.4u±U.ou 


U 


K TT + 7T + 7T~ 


0.60±0.22 


1.20±0.44 





0.60±0.22 


K + 7T + 7T~7T° 


1.05±0.60 


1.05±0.60 


1.05±0.60 


1.05±0.60 


K°7T + 7T°7T 


0.45±0.17 


0.45±0.17 


0.90±0.34 





K+3tt° 


0.15±0.13 





0.45±0.39 





Total 




125.5±11.1 


112.5±8.0 


46.6±6.8 
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Table XVIII: Inclusive yields of kaons from various final states in D s decays. 



Mode 


B{%) 


B(K+)(%) 


B(K°)(%) 


B(K-)(%) 


B{K\%) 


T + U T 


6.6 ±0.6 


0.14+0.01 


0.10+0.01 


0.03 


0.01 


<f>£+v t 


4.72±0.52 


2.32+0.26 


1.60+0.18 


2.32+0.26 


1.60+0.18 


K + K° 


2.98±0.18 


2.98+0.18 








2.98+0.18 


K + K 7T + 


5.50±0.28 


5.50+0.28 





5.50+0.28 





k°k\+ 


5.76±0.96 





5.76+0.96 





5.76+0.96 


k+k\° 


3.56±0.67 


3.56+0.67 








3.56+0.67 




5.6 ±0.5 


5.6 + 0.5 





5.6 + 0.5 





K°K°7i + n 


4.36±1.68 





4.36+1.68 





4.36+1.68 


K + K°7i + n- 


1.92±0.26 


1.92+0.26 








1.92+0.26 


K°K-it + it + 


3.28±0.24 





3.28+0.24 


3.28+0.24 





K + K°n°n 


0.80±0.67 


0.80+0.67 








0.80+0.67 


K + K~7T + 7T + 7T~ 


0.88±0.16 


0.88+0.16 





0.88+0.16 







0.70±0.21 





0.70+0.21 





0.70+0.21 


Ii ./. V /( /I /I 


n 46±0 34 


46+0 34 

W . t:\J _I_W . -J t: 


n 


n 

W 


n 4fi+n 34 


K° K-TT + TT + TT 

11 il /I /I /I 


23±0 16 


n 


n 9S+0 16 


23+0 16 


n 


K+ K~7r + 7T 7r° 

ii ii /I /I /I 


59+0 1 6 


59+0 1 6 


n 


59+0 1 6 


n 


A A 7T 7T 7T 


U.4 i ±U.lo 


U 


U.4 < ±U. io 


U 


n a n i n 1 k 
U.4/ ±U.lo 


A T A n n n 


U.Uo+U.Uo 


U.Uo+U.Uo 


u 


u 


U.Uo+U.Uo 


K°7l + 


0.25+0.03 





0.25+0.03 








K + TT° 


0.08+0.02 


0.08+0.02 











K + T] 


0.141+0.031 


0.14+0.03 











K + i 


0.16+0.05 


0.16+0.05 













0.69+0.05 


0.69+0.05 











K°n + TT° 


0.61+0.35 





0.61+0.35 








K+n°n° 


0.23+0.18 


0.23+0.18 











K°7l + 7Y + 7Y^ 


0.60+0.22 





0.60+0.22 








K + 1T + 1T-1T 


1.05+0.60 


1.05+0.60 











K°n + n°n° 


0.45+0.17 





0.45+0.17 








K+3ir° 


0.15+0.13 


0.15+0.13 











Total 




27.3+1.4 


18.4+2.0 


18.4+0.7 


22.7+ 2.2 
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Table XIX: Inclusive yields of rj, rf, <p, and u from various final states in D s decays. 



Mode 


Y2I07 \ 

B[%) 


Y2l ^\ 107 \ 


Y2l 107 \ 

B{j] )(/o) 


3(0) (%) 


d{uj){7o) 


vrvt 


5.8±1.2 


f O 1 1 o 

5.8±1.2 


A 









A 




1 11+ 


2.04±0.66 


1.33±0.43 


o n a in z^/ -1 

2.04±0.66 







n a /' i n no 

0.06±0.02 


(prv e 


/I 170 i n rn 

4.72±0.52 


n n/" 1 i n ni 

0.06±0.01 


A 




4 


,72±0.52 


A 






4.38±0.35 


n n/" 1 i n ni 

0.06±0.01 





4 


o o i n or 

,38±0.35 


A 




i „+ 
<PP 




n 1 i _i_n no 
0.11±0.03 


n 




8 


43±2.34 


A 






"i on i n oo 

1.20±0.22 


n no 

0.02 





1 


c\r\ \ c\ c\c\ 

.20±0.22 


A 






n on i n 1 r 

0.80±0.15 


n n 1 

0.01 








,80±0.15 


A 






1 r o i n o i 

1.58±0.21 


1 r o t n o 1 

1.58±0.21 


A 









A 






n o^-i-n no 


n 
U 


A 

u 







n o^-i-n no 


1 

VP 


13.0±2.2 


13.0±2.2 













1]'7l + 


3.8±0.4 


2.48±0.26 


3.8±0.4 







0.11±0.01 


r)'p + 


12.2±2.0 


7.97±l.32 


12.2±2.0 







0.37±0.07 


K + r] 


0.141±0.031 


0.14±0.03 













K+rj' 


0.16±0.05 


0.10±0.03 


0.16±0.05 










Total 




32.7±2.9 


18.2±2.1 


19.2±2.4 


0.8±0.1 



For example, both rf{— > 7r + ~tt~ if) p + "(— > 7r + 7r°) and r/7r + cj(^ 7r + 7r~7r°) contain the same 
final-state particles t]2t{ + -k^ -k . In the CLEO paper reporting large branching fractions 
for D s — » r]^p + (where was identified by its two-photon decay mode) [17] . no dis- 
tributions in M(r]) or M(rf) were shown. If a sizable fraction of the measured value 
B(D S — > r]'p + ) = (12.2 ± 2.0)% is, for instance, due to underlying D s — > r]7i + u events, 
then our prediction for inclusive B(r)') becomes smaller while B(u) becomes correspond- 
ingly larger. The effect of such events on B(r]) combines a decrease from a smaller value 
of B(D S — > i]'p + ), where rf decays dominantly to states involving rj (see Table VI), and 
a positive contribution from B(D S — > i]'n- + uj). 

Similarly, the decay mode D s — > ou(—> 7r°7)7r + 7r°, where one photon from 7r° may be 
missing in some events, can mimic D s — > rj(— > r y r y)p + {-^ 7r + 7r ). If a small fraction of 
B(D S — > f]p + ) = (13.0 ± 2.2)% stems from underlying D s — > ujtt + tt° events, then our 
predictions for inclusive B{rf) and B{oj) decrease and increase, respectively, roughly (i.e. 
neglecting differences in efficiencies) by a ratio B(uo — > 7r°7) : B(r] — > 77) = 1 : 4.4. 

• The measurement of B(rj), B(r)'), and B{oj) may help to shed light on specific decay 
mechanisms. For example, the decay D s — > ujti + 7t , represented by the quark annihi- 
lation process cs — > ud, could have a sizable branching ratio. While D s — > ujtt + with 
B = (0.25 ± 0.09)% is forbidden by G-parity and requires preradiation of u or rescat- 
tering [5], D s — > uoti + 'k q can occur through ordinary quark annihilation (as long as some 
mechanism overcomes helicity suppression) , and is expected to have a considerably larger 
decay rate. As a second example, the decay D s — > U7i + r] could arise either from WA or 
from the transition cs — > ss+ (charged weak vector current), where the charged weak 
vector current produces wk + . 

We look forward to experimental tests of these predictions. 
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Note added in proof: Inclusive D s decays have now been reported by the CLEO Collabo- 
ration [19]. They are largely in agreement with the predictions of Tables XVII -XIX, with the 
notable exceptions of B(rj') = (11.7 ± 1.7±0.7)% and B(lu) = (6.1 ± 1.4±0.3)%. Thus, current 
world averages [7] apparently overestimate the sum of D s branching fractions involving rj' and 
greatly underestimate those involving u. 
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